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The phase transition present in the linear-chain and square-lattice cases of the transverse Ising 
model is examined. The extended coupled cluster method (ECCM) can describe both sides of 
' the phase transition with a unified approach. The correlation length and the excitation energy are 

^\ , determined. We demonstrate the ability of the ECCM to use both the weak- and the strong-coupling 

0^ ' starting state in a unified approach for the study of critical behavior. 
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Cn ■ I- INTRODUCTION 

^ • The study of the magnetic properties of materials from microscopic theory has always proved difhcult. Even at the 
' level of such simplified magnetic models as the Heisenberg model, the Ising model, and the Hubbard model, intensive 
' theoretical and numerical work is required to model the interesting counter-intuitive behavior they exhibit. 
\^ I An obstacle to studying quantum magnets, encountered by various theoretical approaches, is due to the uniqueness 
CO i of the state at the phase transition point. Most techniques can describe a particular phase of a quantum magnet but 
' may be unable to describe the other phase, or phases, of the system and, hence, the phase transition itself. Clearly, 
a consequence for such techniques is that conclusive evidence of symmetry breaking cannot be found. 
All techniques dealing with many-body systems experience the difficulty of being able to include sufficient correlations, 
without making prior assumptions about the leading correlations, in order to provide a good description of the 
"j^ important behavior of the system. Unsurprisingly, the subtleties of a quantum magnet at the phase transition, which 
require the inclusion of a great many correlations to model adequately, are notoriously difficult to obtain. Many 
I , theoretical approaches |Q] yield a mean-field like description of the system, which does not reveal the distinctly 
''O ' different behavior present in differing dimensions, contrary to the true physical behavior. |2| 
^ In this paper, we present a study of the transverse Ising model via an improved approach |3,H recently developed for 
^ ' spin systems. This approach involves the application of the quantum many-body technique known as the extended 
coupled-cluster method (ECCM). A priori, we know that the ECCM is potentially a good technique for studying the 
phase transition as it possesses certain features that allow it to study spontaneous symmetry breaking. In particular, 
the ECCM can yield a single solution that correctly describes the physical behavior on either side of the phase 
transition. The more widely applied normal coupled-cluster method (NCCM), which is a restricted version of the 
ECCM, at a given level of approximation, does not possess this feature of the ECCM. A comparison of the results from 
these techniques will demonstrate the superiority of the ECCM with regard to the study of the phase transition. Most 
importantly, we aim to show that the improved approach outlined in Q yields trustworthy and accurate numerical 
data that can provide real insight into the physical behavior of the system. 

The transverse Ising model has two distinct phases; the ferromagnetic phase, and the paramagnetic phase, which 
are tuned by a single coupling constant. On the linear chain, this model has some interesting properties, such as 
the duality, which allows one to relate the ground-state energy at a particular value of the coupling constant to 
the ground-state energy at the reciprocal value of the coupling constant. This property and the fact the model is 
exactly solvable on the linear chain makes it a good test-ground for the ECCM. With the confidence gained from the 
linear-chain model, we will then study the square-lattice case. 
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II. THE TRANSVERSE ISING MODEL 



The transverse Ising Hamiltonian has two competing terms on a lattice composed of spin-half spins; the ferromagnetic 
term, which forces the system to align in the z-direction, and the paramagnetic term, which forces the system to align 
along the positive x-axis: 

dN N 

i/ = -^afa?-A^af , (2.1) 

(iJ> i 

where ct* are the Pauli matrices, A is the external magnetic field, N is the number of lattice sites, d the dimensionality 
of the lattice, and the sum over (i, j) runs over all nearest-neighbor pairs of lattice sites i and j. In the weak-coupling 
limit, A = 0, the spins will align themselves along the z-axis, in either the positive or the negative direction. In the 
strong-coupling limit, A — > oo, the magnetic field forces the spins to align along the positive x-axis. Between these 
two phases a phase transition occurs for some critical coupling Ac. 

Starting from the strong-coupling limit, the Z(2) symmetry is preserved, which results in the (ct^) = Af^ being zero. 
Beyond the critical point this symmetry is broken, and two degenerate ground states exist, with ±M^. Any ab-initio 
method would have difficulties in breaking such a symmetry. We will show that within a specific approximation, we 
can start with an uncorrelated model state in both limits, which can either incorporate the Z{2) symmetry or break 
it, in the latter case two exactly degenerate states exist. 



A. Exact results on the linear chain 



On the linear chain the model is exactly solvable. [^,^ After a Jordan- Wigner transformation the Hamiltonian can 
be diagonalized using a Bogoliubov transformation. The ground-state energy is given by: 



- / dq^l + 2Acosg A^ 

Jo 



N 



-3(1 + ^)^(7^^) ' (2-2) 



where f is a complete elliptic integral of the second kind. This energy satisfies the interesting duality, as could be 
concluded from the duality of the original Hamiltonian: 



Eg{X) = XEg{X-') . (2.3) 



The magnetization M°' — (ct") is given by: 



Jo V 1 + 2A cos q + X'^ 
= 6(1 - X) {I - X^)^ , (2.5) 

where 9{x) is the step function; 9{x < 0) = and 9(x > 0) = 1. The transverse magnetization can be related to 
the derivative of the ground-state energy, however, the longitudinal magnetization is derived from the limit of the 
spin-spin correlation function with a lattice size taken to infinity, therefore the vanishing long-range order for A > 1 
is a subtle result. |^ 

B. The mean field 

A parameterization of the wave function is given by: 

f-sme/2\ f~sm9/2\ (~sm9/2\ f-sin9/2\ \ 
^ \ cos9/2j\ cos9/2j'y cos9/2j\ cos9/2j'"/ ' ^ °'' 

From the variational principle the mean-field solution yields: 
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e = -arcsinA ; A < 1 , (2.7) 

= -| ; A>1 , (2.8) 

^ = -d(l + P) ; A<1 , (2.9) 

= -2dA ; A > 1 , (2.10) 

M'' = 6* (1 - A) Vl - A2 , (2.11) 



where A = 2d\. Note that 9 — tt — 9 is also a solution with the same energy but opposite magnetization. This is 
the Z{2) symmetry partner of the ferromagnetic ground state. 



III. THE EXTENDED-COUPLED CLUSTER METHOD 



The Coupled-Cluster Method (CCM) which exists in two distinct versions, the NCCM and the ECCM, is a powerful 
technique for quantum many-body calculations. Only the essential features of the CCM will be reviewed here, further 
details can be found in the literature. @,^,D There are four key features of the CCM. Firstly, the model state |$), a 
uncorrelated state, which is usually the true ground state in a certain limit, serves as the starting state. Secondly, 
the correlations are incorporated in the the ket state in the exponentiated form: 

|vl/) = e^|<i>) , (3.1) 

where the correlation operator S is composed solely of a complete set of mutually commuting multiconfigurational 
creation operators that annihilate the bra model state ($|. The exponential form ensures the correct summation of the 
independent correlations in the calculation, and spreads correlations over the whole system although the correlations 
in S might be only local. 

The third key feature of the CCM appears in the parameterization of the bra state {^\, 

{^{1 + S)e-^ ^ {^e^" e-^ , (3.2) 

where the first parameterization gives rise to the NCCM, and the second parameterization gives rise to the ECCM. 
The importance of both of the above bra-state parameterizations are that they allow the similarity transform to 
be incorporated into the CCM calculation of any observable. Primarily, the similarity transform only yields terms 
for which the correlations present in S are linked to the observable, and, hence, subsequent expectation values are 
size-extensive quantities. Importantly, the presence of the similarity transform means that no further truncations are 
necessary than the truncations of the number of configuration of creation operators in S. In this way the similarity 
transform can be seen as the defining feature of the CCM, since its absence would render an approximat e ca lculation 



impractical. For example, if one would attempt a variational calculation with the ket state of Eq. (BT) and its 
Hermitian conjugate, several further approximations are required to calculate the functional. Therefore, the actual 
ket state obtained from the same parameterized variational wave functional depends on the associated bra state and 
the corresponding variational principle. 

This brings us to the fourth key feature of the CCM, which is related to the parameterization of the bra state in Eq. 

( |3.2D and concerns the terms arising from e'^ and S. These terms are present in the ECCM and NCCM respectively 
and are solely composed of destruction operators; the Hermitian conjugates of the creation operators in S. Therefore, 
the functional is a polynomial in the bra and ket coefficients, associated with the strength of the different 

creation and annihilation operators in S and S respectively. 

Both of the NCCM and ECCM parameterizations result in the maintenance of some useful formal properties, which 
remain valid in any approximation scheme, e.g., Feynman-Hellman theorem and the linked-cluster theorem for the ket 
state coefficients. ||] However, the ECCM, unlike the NCCM, parameterization retains, in any approximation scheme, 
the full cluster separability: 

lim (ArBr') - (^r)(Sr') , (3.3) 
|r — r I — >oo 

where Ay. and B^i are localized operators acting at positions r and r , respectively. Therefore, it allows us to separate 
the correlation spin-spin function {cr^cr^,) into a size-extensive part, the long-range order; (cr^)^ = (M^)^, and the 
true correlation function that vanishes at infinite distance. 
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A. Transformation of the Hamiltonian in terms of a canted model state 



For a unified description of the different model states |<I>), we perforin a unitary transformation such that the model 
state, which is a mean field state where the spins point in a particular direction, is a state composed of down-pointing 
spins. The details of this transformation have been given elsewhere. Hence, the creation operators in S will be 
products of cr"*"'s at different sites. Using the rotation matrix U, for this purpose, defined by, 

U = exp (^i9^^ = cos 1 + * sin (J^ , (3.4) 

one can obtain the transverse Ising Hamiltonian defined with respect to an arbitrary canted model state, in terms of 
rotated operators corresponding to the down-spin model state. An illustration of the particular canted states of the 
ferromagnetic and paramagnetic states, described by the angle 6, are respectively given by: 

= ; Uiii---)A=o , 9 = -n/2 ; | ^^^^ • • •)a^oo • (3.5) 
The rotated transverse Ising Hamiltonian is of the form: 

dN 

UHW = - ^ [sin^ Oaf erf + cos^ Oa^cr? + sin 20crfcr?] 

(iJ> 

N 

- A ^ [- sin 9(7? + cos Oa?] . (3.6) 

i 

In most cases we use only two rotations: the ferromagnetic Hamiltonian, with 9 — 0, and the paramagnetic Hamilto- 
nian, with 9 = 



B. The SUBl scheme and symmetry breaking 



The ECCM SUBl approximation scheme only retains the one-body correlations present in 5* and S in Eqs. (3.1) 



and (3.2) defined with respect to the rotated model state, given by: 

S^kJ2a+ ; S"^k"j2ar . (3.7) 

i i 

Performing the CCM SUBl scheme is equivalent to a mean-field calculation with an important proviso. Crucially, 
the ECCM SUBl scheme can yield the full mean-field ground-state solution, however, the NCCM SUBl scheme can 
not yield the ground-state solution for the region beyond the phase transition. This feature of the NCCM has dire 
implications for the study of symmetry breaking. It is a little known fact that, in practice, the CCM yields a solution 
that preserves the shared symmetry of the model state and the Hamiltonian Fortunately, via the use of the ECCM 
SUBl scheme, it is possible to obtain a symmetry-broken solution in a correlated ECCM calculation. 
The relationship between the mean-field ground state and the ECCM one-body coefficients defined with respect to 
the canted model state, where the spin are rotated downwards, is given by: Q 

fc = -tan^0 ; k" = -^sm9 . (3.8) 

In a correlated ECCM calculation, these relations are only valid in the classical regimes, A = and A = cxd, of the 
system. By using the values of k and k , from Eq. ( |3.8D , starting from an initial point in the classical region of 
the phase where the chosen model state is not the classical wave function, one can determine the symmetry-broken 
solution. 

The above ansatz is put into practice to study the possibility of broken Z(2) symmetry in the transverse Ising model. 
Employing the paramagnetic model state, with the values k = 1 and k — 1/2, which has the advantage of not 
artificially breaking the Z{2) symmetry, corresponds to the ferromagnetic state. At the ferromagnetic point A = it 
yields the symmetry-broken solution. For k — and k =0 the symmetry preserved solution can be determined, for 
which Af^ = 0, from the limit X ^ oo. 
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C. The diagrammatic representation of the SUB2-n scheme functional 



In order to perform the numerical ECCM calculation efficiently, a diagrammatic representation is implemented. 
The ECCM SUB2-n approximation scheme is implemented here and is described in greater depth in References 0,1). 
Essentially, the SUB2-n scheme retains only one and two-body correlations, where the truncated correlation operators 
are of the form: 

5 = /C ^ (7;+ + ^ 6x(r)fT+fTi+r , S" = k' ^(^i K{v)^i ''i+r ' (3-9) 

i i,r i i,r 

where x(r) yields a label for each vector distinct under lattice symmetries. 

The representation of the combinations of operators which give rise to various physical quantities is outlined in 
Reference . The transverse Ising functional, defined with respect to the ferromagnetic model state arising from the 
ECCM SUB2-n scheme is diagrammatically expressed in Fig. |l|. 



D. The correlation length and the excitation energies 



Although the order of the approximation schemes is not sufhcient to study the long-distance behavior of the correlation 
function, we can examine the correlation length: 

^ = ^(E-f-r)-]^(E-f)' ■ (3-10) 

A similar correlation length for the z-components of the spins is computationally too involved to be of any practical 
use. This expression for the ^ is easily motivated in the linear-chain case, where 

((Tf (jf ) OC L-l'-jl/'^ + {A-Pf (3.11) 
where is the long-range order, such that ^ = (e~ — 1)^^ ~ C- 

The excitation energies are calculated using the random phase approximation (RPA), which means that we assume 
that the excitations are small, harmonic fluctuations around the ground-state solution. This leads to a linear eigenvalue 
problem, which is solved numerically. We can study the RPA spectrum, since the second-order derivatives are 
available for the functional, which is stored quasi-analytical for the high-order numerical calculations. At the phase 
transition the excitation energy goes to zero for a second-order phase transition, however, we know from NCCM |^ 
that the vanishing excitation energy is also the generic behavior at a termination point. Similar behavior could occur 
for ECCM. 



IV. NUMERICAL RESULTS 



Usually, there is not much to be gained from studying the ground-state energy in great detail with the purpose of 
examining the performance of a method. Its convergence is much faster than that of any other quantity, and generally 
it is not an indication of the accuracy of the ground-state wave function. Moreover, for the transverse Ising model 
there is no value A, where the numerical value of the energy is of specific interest. The position of the phase transition, 
Ac, the only candidate for a point of interest, is not exactly known. The interested reader can contact us for precise 
numerical values for the ground-state energy. Instead, we focus on the excitation energy and the magnetization. 
It is important to note that there are four specific cases studied here. Two different model states: the symmetry- 
broken, down-pointing ferromagnetic model state and the symmetry-preserved, sideways-pointing paramagnetic model 
state. These are two distinct approaches to the problem, which lead to two different numerical calculations. However, 
both can be studied, starting from both of the classical limits, A = and X — oo. The one-body terms, from the SUBl 
calculation, allow us to rotate each state from the model state to the starting state, which leads to some unexpected 
differences, which are summarized in Table |. 

The ferromagnetic model state, which breaks the Z{2) symmetry, generally does much worse that the paramagnetic 
model state. We consider the two classical states that can arise with the paramagnetic model state. Firstly, the 
starting state at A = breaks the Z{2) symmetry since it has = 1. However, following this solution for the 
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ground state for increasing A, it turns back and connects up with the symmetry partner; —M^, and ends at 

the starting value A = 0, but now at with = — 1. At low orders this is a smooth transition, where goes to 
zero, however, at high orders it jumps at a finite value of to —M^. At this point the excitation energy vanishes, 
for any order of approximation. Therefore both symmetry partners of the ground state are exactly degenerate and 
part of one continuous solution. 

Secondly, starting with the paramagnetic model state at large values of A, wc find that the remains zero, so the 
symmetry is never broken. This solution crosses the symmetry-broken solution as the latter turns back at Af^ — 
and Ac. However, looking at the excitation energy of the symmetry preserving solution, it tends to zero on the square 
lattice (see Fig. H), but remains finite on the linear chain (see Fig. ||), signaling different behavior. Moreover, the 
ground-state solution continues beyond the point where the excitation energy goes to zero, therefore this behavior 
is not associated with the possible generic behavior at the termination point. The value of Ac, where the excitation 
energy vanishes, is lower than the termination point of the solution with the paramagnetic model state starting from 
A = 0. 

The results from the square-lattice case (see Fig. ^ and Tabl e O ) are in greater agreement with those from other 
methods than the linear-chain results (see Fig. ^ and Table are with the exact results. Although this is to 
be expected, as the large quantum fluctuations on the linear chain is numerically difficult to deal with. Our most 
accurate means of determining the critical coupling constant is with the paramagnetic model state starting from the 
ferromagnetic limit A = 0. It yields a value of Ac = 3.014 which is close to the best high-temperature series expansion 
results of Ac = 3.0441(4) and the best low-temperature series expansion results ||ll[ of Ac — 3.041(3). Note that 
the series expansions, respectively, use the ferromagnetic and the paramagnetic model state. Other methods find 
similar, although less accurate results. |l2| The details are summarized in Table ||. 

To investigate the model-state dependence further, we look at angles 9 in between the ferromagnetic angle 6 — and 
the paramagnetic angle 9 = The largest angle at which the does not jump is rather close to the ferromagnetic 
angle. The for these different model states are shown in Fig. ^ Note that this solution continues to A ^ oo 
because the Z{2) symmetry of the model state is broken. A continuous solution that traverses across the phase 
transition might be in high demand, however, the solutions, from a canted angle model state, do not have the right 
properties as the excitation energy becomes complex before the phase-transition point. 

The correlation length increases in the vicinity of the phase transition (see Fig. and shows good convergence. 
However, whether the correlation length is diverging at the phase transition is not clear. For the linear-chain case 
there are some indications of a sudden increase. More interestingly, for lower-order approximations, the maximum 
correlation length occurs before the termination point, and, as the order increases, the termination point and the 
extremum in the correlation length converge. 



V. CONCLUSIONS 



The purpose of studying the transverse Ising model with the ECCM was to gain further knowledge of the optimal 
approach for its implementation before tackling more complicated physical models. Two major considerations came 
to light: the pivotal importance of the model state in the CCM calculation and that, importantly, the ECCM is 
the version of the CCM which can show symmetry breaking. Encouragingly, the ECCM SUB2 numerical results are 
accurate in comparison to other numerical techniques and is in agreement with the series expansion result [ p^|Jll|| for 
Ac on the square lattice. 

The calculations in this paper necessitates the inclusion of the full ECCM SUBl part of the functional, something 
we avoided in earlier work Q| because of the computational difficulties. However, the advantages of this inclusion are 
clear, since it allows us to study both phases in one consistent approach. Unfortunately, the order of approximation 
to which this calculation can be pursued is lower, because of the greater complexity. However, in the case of the 
transverse Ising model high-order approximations are not important since the convergence is good, and apart from 
the results for correlation length we have fully converged results within the SUB2 approximation, at SUB2-20. 
The fact that the optimal ECCM SUB2 solution terminates in the critical region reveals that further improvements 
are necessary. Using more complicated model states, i.e., a dimerized paramagnetic model state, could improve the 
situation. It seems that ECCM can only be implemented efficiently in the SUB2 truncation scheme, which is a serious 
drawback. However, a new approximation scheme has recently been developed, [|| which employs "block-spin model 
states" within the auspices of the SUB2 scheme. Advantageously, this scheme can be viewed as a way of performing 
a re-summation of the diagrams arising from the usual single-spin model states and as such can allow the inclusion 
of a great-many correlations. The formulation of an algorithm |^ for this scheme further increases its viability via 
computerization. 
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FIG. 1. The diagrams arising from the transverse Ising model functional with the ECCM SUB2-n scheme, which yields the 
ground-state energy, where d is the dimensionality of the hyper-cubic lattice and A is the external magnetic field. 
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FIG. 2. The excitation energy, _Ea;, as a function of the transverse magnetic field, A, on the square-lattice case from the 
ECCM SUB2-n scheme with the paramagnetic model state and a mean-field calculation. 
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FIG. 3. The excitation energy, Ex, as a, function of the transverse magnetic field. A, on the linear-chain case from the ECCM 
SUB2-n scheme with the paramagnetic model state and a mean-field calculation. 
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FIG. 4. The longitudinal magnetization, M^, as a function of the transverse magnetic field. A, on the for the linear-chain 
case from the ECCM SUB2-n scheme with the paramagnetic model state, the highest-order NCCM LSUBn scheme with 
the ferromagnetic model state, a mean-field calculation, and the exact solution. The LSUBn truncation scheme retains only 
multiconfigurational creation operators in a localized area of n contiguous sites. 
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FIG. 5. The longitudinal magnetization, M^, as a function of the transverse magnetic field, A, on the for the linear-chain 
case from the ECCM SUB2-n scheme with the paramagnetic model state, the highest-order NCCM LSUBn scheme with the 
ferromagnetic model state, a mean-field calculation, and the exact solution. 
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FIG. 6. The longitudinal magnetization, Af^. as a function of the transverse magnetic field. A, on the square-lattice case 
from the ECCM SUB2-n scheme with the paramagnetic model state, a critical canted model state, which allows the solution 
to be continuous, and the ferromagnetic model state. 
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FIG. 7. The correlation length ^ as function of the coupling constant A, for both the linear-chain case (left) and the 
square-lattice case (right) with increasing orders of approximation. For the linear-chain case the correlation length is converged, 
however, not completely so for the square-lattice case. For the square-lattice case, each line represents an approximation, where 
all correlators in a square, with lengths: 1, 2, 3, and 4, are taken into account, corresponding to 2, 5, 9, and 14 distinct 
correlators. 



TABLE I. A summary of the main features for the square-lattice case, with different model states and different starting 
points. The qualitative difference with the linear-chain case is only for the excitation energy for the paramagnetic model state 
starting from large A. 



starting value 


paramagnetic model state 


ferromagnetic model state 


A = 


solution terminates 


solution continues 


AT 


Ex vanishes 


Ex turns complex 




terminates 


decreases 


A = oo 


solution continues 


solution continues 


Ai 


Ex vanishes 


Ex turns complex 




= 


increases 



TABLE IL A comparison of the critical value of the external magnetization at which various techniques determine the 
presence of a phase transition on the square lattice. 





VQMC^ 


HNC CBF^ 


HT Series Expansion'' 


LT Series Expansion"^ 


ECCM SUB2 


Ac 


3.15 ±0.05 


3.12 


3.0441(4) 


3.041(3) 


3.014 



^ Reference 
^ Reference |TO 
Reference [O 



TABLE III. A comparison of the critical value of the external magnetization at which various techniques determine the 
presence of a phase transition on the linear chain. 
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Exact 


ECCM SUB2 


Ac 


1.206 


1.22 
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1.12 



Reference fl^ 
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